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Abstract 

Let P = circ(Pi,P2, . . . , P n ) and Q = circ(Qi,Q2, • • • , Qn) benxn 
circulant matrices where P n and Q n are nth Pell and Pell-Lucas numbers, 
respectively. The determinants of the matrices P and Q were expressed 
by the Pell and Pell-Lucas numbers. After, we prove that the matrices P 
and Q are the invertible for n > 3 and then the inverses of the matrices P 
and Q are derived. 



1 Introduction 



Circulant matrices have a widely range of application in signal processing, cod- 
ing theory, image processing, digital image disposal, design of self-regress and 
so on. Also, numerical solutions of the certain types of elliptic and parabolic 
partial differential equations with periodic boundary conditions often involve 
linear systems Cx — b with C a circulant matrix [9-11]. 

The n x n circulant matrix C' n = circ(co, c±, . . . , c„_i), associated with the 
numbers cq, c\, . . . , Cn-i, is defined by 



C Cl 
Cn-l C 



C-2 
Cl 



C3 



C n -2 C„_i 
Cn-3 C„_2 



C() 
Cn-l 



Cl 
C() 



(1) 



Circulant matrices have a wide range of applications, for examples in sig- 
nal processing, coding theory, image processing, digital image disposal, self- 
regress design and so on. Numerical solutions of the certain types of elliptic 
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and parabolic partial differential equations with periodic boundary conditions 
often involve linear systems associated with circulant matrices [9,10]. 

The eigenvalues and eigenvectors of C n are well-known [4,14]: 

n-l 

Aj = ^CfeW 5fe , j = 0, l,...,n- 1, 

fc=0 

where cj := cxp(^p) and i := \/— T and the corresponding eigenvectors 

X j = (l,LJ j ,CJ 2j ,...,J n - 1 »), J =0,1,. ..,71-1. 

Thus we have the determinants and inverses of nonsingular circulant matrices 
[1,3,4,14]: 

n—1 n—1 

det(C„) = II(E c ^ fe ) 

j=0 fc=0 

and 

C" 1 = drc(a ,ai, . . . ,a„_i) 

where ak — ^ Y^j=o Aj w ~" ,fc , and fc = 0, 1, . . . , n — 1) [4]. When n is getting 
large, the above the determinant and inverse formulas are not very handy to use. 
If there is some structure among Co, c±, . . . , c„_i, we may be able to get more ex- 
plicit forms of the eigenvalues , determinants and inverses of C n . Recently, stud- 
ies on the circulant matrices involving interesting number sequences appeared. 
In [1] the determinants and inverses of the 
circulant matrices A n = circ{F\, F2, . . . , F n ) and B n — circ{L\, L2, ■ ■ ■ , L n ) 
are derived where F n and L n are nth Fibonacci and Lucas numbers, respec- 
tively. In [2] the r-circulant matrix is defined and its norms was computed. 
The norms of Toeplitz matrices [??] involving Fibonacci and Lucas numbers are 
obtained [5]. Miladinovic and Stanimirovic [6] gave an explicit formula of the 
Moore-Penrose inverse of singular generalized Fibonacci matrix. Lee and et al. 
found the factorizations and eigenvalues of Fibonacci and symmetric Fibonacci 
matrices [7]. 

The Pell and Pell-Lucas sequences {P n } and {Q n } are defined by P n = 
2P„_i + P n -2 and Q n = 2Q„_i + Q„_2 with initial conditions P = and 
Pi = 1 and Q = 2 and Qi — 2 for n > 2, respectively. The Pell numbers are 
used to approximate the square root of two, to find square triangular numbers, 
to construct integer approximations to the right isosceles triangle, and to solve 
certain combinatorial enumeration problems. Let P = drc(Pi, P2, . ■ . , P n ) and 
Q = circ(Qi,Q2, • ■ • , Q n )- The aim of this paper is to establish some useful 
formulas for the determinants and inverses of P and Q using the nice properties 
of the Pell and Pell-Lucas numbers. Matrix decompositions are derived for P 
and Q in order to obtain the results. 
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2 Determinants of circulant matrices with the 
Pell and Pell-Lucas numbers 



Recall that P = circ(Pi, P2, . . . , P n ) and Q = circ(Qo, Q\, . . . , Q n -i), i.e. where 
Pk and Qfc are Pell and Pell-Lucas numbers, respectively, with the recurence 
relations P& = 2Pfc_i +Pk-2, Qk = ^Qk-i+Qk-2, the initial conditions P = 0, 
Pi = 1, Q = 2 and Q x = 2 (fc > 2). Let a and /3 be the roots of x 2 - 2x- 1 = 0. 
Using the Binet formulas [11, p. 142] for the sequences {P n } and {Q n }, one has 

and 

Q„ = (3) 



Theorem 1 Let n > 3. T/ien 



n-l 



det(P) = (P - P i+1 )"- 2 (P! - 2P„) + J2 [Pk-iP,:- k (Pi - Pn+i f- 2 ] (4) 
where Pk is kth Pell number. 



k=2 



Proof. Obviously, det(P) = 104 for n = 3. It satisfies ((H). For n > 3, we selec 
the matrices M and N so that when we multiply P with M on the left and TV 
on the right we obtain a special Hessenberg matrix that have nonzero entries 
only on first two rows, main diagonal and subdiagonal: 



and 



M := 



N := 



1 

-2 

-1 





1-2 
1-2-1 



Pn 


Pi 


~Pr 
P n 


+ 1 


Pi 


-P, 
P n 


+ 1 



Pl-PnA 



n—4 



(pi-k+i) 
1 





1 

1 
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Notice that we obtain the following equivalence: 



S 



MFN 



f/n 


Pfi-1 


Pn-2 


Pn-3 


Pi 


9n 


Pn-2 


Pn-3 


Pn-A 


Pi 




Pi - Pn+1 










-Pn 


Pi - Pn+1 











-Pn 


Pi — Pn+1 





-Pn Pi — Pn+1 

— Pn Pi — Pn+1 



and S is Hessenberg matrix, where 

n-l 



Pr, 



g' n := Pn + J2 Pk (i 

k=2 ^ 
n-l 

g n :=P 1 -2P n + Y / Pk-i . p p 

. * i ±i 



— Pn+1 
Pn 



k=2 



n+1 



n—k 



Then we have 

det(S) = det(M) det(P n ) det(iV) = (P 1 - P n+1 ) n - 2 g„ 

Since 



for all n > 3, 
and ((IJ follows 



det(M) = det(AT) = 

det(M) det(iV) = 1 



1 , n — 1 or 2 mod 4 
1 , ri = or 3 mod 4 



Theorem 2 Let n > 3. T/ien 



= 2(2-Q„ +1 )"- 2 (2-3Q„)+2 £ [(Q fe+1 - 3Q fe )(Q„ - 2)"^ fc (2 - Q n+1 ) fc - 

fe=2 

(6) 

where Qk is kth Pell-Lucas number. 

Proof. Since n > 3, det(Q) = 2464 for n = 3, satisfies ©■ For n > 3, we select 
the matrices if and L so that when we multiply Q with K on the left and M 
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on the right we obtain a special Hessenberg matrix that have nonzero entries 
only on the first two rows, main diagonal and subdiagonal: 



(7) 



and 



We have 

U = 



Qi 





1 








.. 










-3 








.. 





1 




-1 








.. 


1 


-2 


K := 











.. 


. -2 


-1 










1 


-2 .. 













1 


-2 


-1 .. 









L := 





-2+Q„ 

Ql-Qn + l 

-2+Q n 

Ql-Qn + l 



n-3 



( -2+0, y 

\Q1-Q„ +1 J 

( -2+Q„ \ 
\Ql-Q„+i J 

1 



u 'n Qn-1 Qn-2 

M n Qn 3Q n _i Qn—1 3Qn — 2 
Ql — Qn+l 

2 — Qn Ql — Qn+l 

2-Qn 



and U is Hessenberg matrix, where 





1 

1 








Qz Q 2 
Qi — 3Q3 Q3 — 3Q2 



Ql — Qn+l 

2 — Qn Ql — Qn+l 



■■= Qi - 3Q„ + (Gfc+i - ( Ql Q " Qn 2 +1 ) 



Then we obtain 

det(I7) = det(A-) det(Q) det(L) = Qi(Qi - Q n +i) n ~ 2 u, 
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Since 

a ±(ts\ i ut\ f 1, n = 1 or 2 mod 4 
det(tf) = det(L) = | _ 1; n ^ 0or3 mQd4 

for all n > 3, 

det(if) det(X) = 1 

and we have ©. ■ 

3 Inverses of P and O 



We will use the well-known fact that the inverse of a nonsingular circulant matrix 
is also circulant [14, p.84] [12, p.33] [4, p.90-91]. 



Theorem 3 Let P = circ(P\, P2, ■ ■ ■ , P n ) is {-avertible when n > 3. 



Proof. From Theorem 1, we have det(P) = 104 ^ and det(P) = -18560 ^ 
for n — 3 and n — 4, respectively. Then P is the invertible for n — 3, 4. Let 
n > 5. TheBinet formula for Pell numbers gives P n — - — 

2, a/3 



u(uj k ) 



-1 and a- (3 = 2^2. Then we have 

'a r -P r 



2 J^ where a + (3 — 



r—k 



r=l 
1 

1 

2y/2 
1 - Pn+1 



E 

r=l 



kr — k 



2V2 

a{l-a n ) ,8(1 -0") 



n 



kr — k 



1 — auj k 1 — f3uj k 



2V2 • , 

v r— 1 

, (l-aw fe ,l-/3w fc 7^0) 



(a - /3) - (a r 



/3" +1 )+a/3w' £ (a™-/3") 



P„w A 



2fc 



1 — auj k — f3u k + ctf3uj 
: = l,2,...,n- 1. 



2/r 



1 - 2w fe - cj 

If there existed w fc (fc = 1, 2, . . . , n — 1) such that u(u k ) = 0, then we would 
have 1 - P n +i - P n oJ k = for 1 - 2u k - uj 2k ^ 0. Hence ui k = 1- p^ +1 ■ It is well 
known that 



exp 



2fc7ri 



[2k^\ [2kTT 



(8) 



where i := v^l. Since 



1--Pr. 



-1 for < 



•2kir 



is a real number, sin (^p-) = so that 



< 2tt. However x 



T is not a root of the equation 



1 — P n +i — P n x — (n > 5), a contradiction, i.e. u(ui k ) 7^ for any w fe where 
fc = 1, 2, . . . , n — 1 and n > 5. Thus, the proof is completed by [1, Lemma 1.1]. 
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Lemma 4 Let C = (cy) be an (n — 2) x (n — 2) matrix defined by 



Pi - P n+ i,i = j 

c ij = { ~ P ni i=j + l 

0, otherwise. 



Then C 1 = (q) is given fry 



c' 



p r 3 • ^ ■ 
(Pi-p„+i)'-^'* - •? 

0, otherwise. 



Proof. Let ^4 := (a^) = CC x . Clearly, = X)fc=i c ik c 'kj- When i — j, we 
have ^ 

a-ii = (Pi - Pn+l)- n 5 = 1- 

If i > j, then 

pi— 3 — 1 pi— 3 

(JWn + i) ,„ " =0; 



"(Pl-Pn+l)*^ ^ (Pi "Pn+l)*" - 

similar for i < j. Thus, CC -1 = I„-2- ■ 



Theorem 5 Let i/ie matrix P 6e P = circ(P\, P2, ■ ■ ■ , P n ) (n > 3j. TTien the 
inverse of the matrix P is 



where 



Pi = 



- 1 = circ(p 1 ,p 2 , . . . ,p n ) 



1 / , 2P„"- J , ^ P n-kP«- 



i + f£» + Y 



9n \ (Pi - Pn+l) n - 2 (Pi - Pn+l)* 



1 / n ~ 2 P pfc-1 

1 / 01 \ " r n-):-l r ii 

pi-3 



■„(Pl - P„+l)*- 2 

n—k 



, i = 3,4,...,n 



/or g n = Pi- 2P„ + P fc -i (p-^) 
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Proof. Let 



U = 



~9n 
1 






3n 
9n 








Pn-2 — Pn-1 


^Pn-3 — Pn-2 
P„-3 


. §*P1-P 2 


g n 
1 





Sr. 







1 




















1 



n—k 



and i? = diag(l,g n ) where ^ = P„ + ^£=2 Pfe ^_p" + J and .g„ = P 1 

n-l ( P \n—k 

2Pn + J2k=2 P k-1 { p 1 -P n+1 ) ■ Then We Can Write 

M¥NU = H®C 

where H C is the direct sum of the matrices H and C. Let T = NU. Then 
we have 

P- 1 =T(H- 1 ®C~ 1 )M. 

Since the matrix P is circulant, its inverse is circulant from Lemma 1.1 [1, 
p.9791]. Let 

P _1 = circ(pi 7 p 2 , ■ . . ,p n ). 
Since the last row of the matrix T is 

Pn-2 Pn-3 Pl-4 



0,1,- 

9n Qn Qn 

the last row components of the matrix P _1 are 

-2 

t / \ 

P2 



P 2 Pi 



•jn yn 



l_ I Pn-k-lPt 1 
9n{ + fe(A-P„ +1 )\ 



P3 
P5 



1 



ffn(Pl — Pn+l) 
P» 

"fln(Pl -P i +l) i 



- E 



f« V^i (Pi " p "+ i)fe 

n-2 



*E 



p 3 - fe p 



fc-i 



Pi 



£ (ft - W 



Pl — Pn+l 



1 P n _fc_lft 



fe-1 



n— 3 



Pl ~ 9n[ + (Pl-P„+l)™- 2 



n — 3 

^E 



Pn-fe- 2 P„ fe - 1 



ftv 



^+l) fe 
-feP^ 1 



E 



Pn-k-^Pn 1 



^(ft-ft^ 
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where g n =P l - 2P n + Y2=l ft-i (n^l)" " ■ If ^ = ElU fp7^> 
(r = 1, 2, . . . , n — 2), then we obtain 



C(2)_ 2 c(l) l '" 



(Pi - P n +l) 2 

and 

pr+l 

g(r+2) oo(r+l) g(r) _ J n _ i n „_4 

i r 1 ■'n+lj 

Hence, we have 

P- 1 = — arc ( 1 + 2S, ( ™- 2) + S^ 3 ) , -2 + S^ 2 ) , -S« , -5,( 2) + 2S« , 
.'/<• v 

o(3) _ 00(2) _ o(l) 9(n-2) _ oo(n-3) _ o(n-4)\ 



1 Clrc ( 1 + 2P r 3 y? Pn-^- 1 _ 2 ^ Pn- k -!P, k - ' 

-arc ii+ _ p v n _ a +2^ rR - P_.^*' + 2^ 



3« V (Pi - P„+i)"- 2 £j (Pi - P„+i) fc S ( p i - 

i p„ pI Pir 3 



k ' 



\n-2 



(P!-P n+1 )' (P!-P n+1 ) 2 ' (P 1 -P (1 + 1 ) 3 '-"' (Pl-P n+1 ) 

by the nice properties of Pell numbers. ■ 

A Hankel matrix A = (ay) is an n x n matrix such that a^j — dj-xj+i. It 
is closely related to the Toeplitz matrix in the sense that a Hankel matrix is an 
upside-down Toeplitz matrix. 



Corollary 6 Let the matrix M be as in f5|). Then its 



inverse is 



Ar 1 = 



1 

B C 



where the matrices 

B=(P n P n -! P„- 2 ... P 2 ) (n _ x) xl 

and C is an (n— l)x (n— 1) Hankel matrix that it has the first row [P n _i, P n -2, ■ ■ ■ > Pi 
and i/ie Zasi column [Pi, 0, . . . , 0] T . 

Proof. The matrix M _1 is obtained easily by applying elementary row 
operations to the augmented matrix [M : I n ]. ■ 

Theorem 7 The matrix Q — circ(Qi, Q2, ■ ■ ■ , Qn) is invertible when n > 3. 
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Proof. Wc show that det(Q) = 2464 f and det(Q) = -1247232 f by 
Theorem 2 for n = 3,4. Then the matrix Q is the invertible matrix for n = 3,4. 
Let n > 5. The Binet formula for Jacobsthal-Lucas numbers yields Q n = a n + 
f3 n , where a + (3 = 2 and a/3 = — 1. Then we have 



v{uj k ) = Q r u kr ~ k = {a r + (3 r ) u k 

(l-aw*,l-0w* t^O) 



r— 1 r— 1 



a(l-a") , P{l-n f , „, , , . , 



1 — aw fc 1 — /3cj fe 
(a + 0) - (a n+1 + (3 n+1 ) + a(3uj k {a n + /T) - 2a0w* 



1 - (a + (3)uj k + a/3uj 2k 

2 - Qn+l + (2 - Q n )0J k 



l-2w fc 



fc = 1,2,. ..,ra- 1. 



We are going to show that there is no uj k , k — 1, 2, . . . , n— 1 such that w(oj fc ) = 0. 
If 2 - Q n+ i + (2 - Q n )uj k = for 1 - 2uj k - uj 2k ^ 0, then w fc = %^ would 
be a real number. By ([5]) we would have sin (^-) = so that ui k = —1 for < 
2^ < 2ir . However x = — 1 is not a root of the equation 1 — Q n+ i + (2 — Q n )x = 
(n > 5), a contradiction, i.e. v{uj ) ^ for any uj k , where k — 1,2,..., n— 1 and 
n > 5. Thus, the proof is completed by [1, Lemma 1.1]. ■ 

Lemma 8 If the matrix A — (dij)"7—i is of the form 



3l - Qn+l, i = j 

aij = { 2 - Q n , i = j + l 
0, otherwise, 



then A 1 = (a^)™-^ is given &y 



{ (Q n -2)*-i 
a[ 3 = <| (Qi-Q.+i)-^ 1 



0, otherwise. 
Proof. Let £? := (6^) = AA^ 1 so that 6.y = 53fc=i a ik a kj- Clearly 
bu = (Qi — Qn+i)- 



Ql — Qn+l 

If i > j, then 

bij = y^fr— 1 0<ik a kj = a i,i—l a i—l,j + a ii a ij 

= (2 _ Qn) (Qi - qIi)^ + {Ql ~ Qn+l) (Qi Q -Qn+iy-^ = ° ; 



similar for i < j. Thus A = /„_ 



2- 
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Theorem 9 Let n > 3. The inverse of the matrix Q is 

(Qr 1 = circ(q 1 ,q 2 ,...,q n ) 

where 



1 L 8(Q n -2) r 
qi — — ii — 



(Ql - Qn+l) n - 



+ E 



(Qri-fc+2 — 3Qn-k+l)(Qn — 2) fc 1 



fc=l 



(Ql-Qn+l) fe 



n-2 



<?2 = —(s + ^Qn-k+l-SQn-k)^ 2)k " 
Un \ fc=l 

4(Q„ - 2) m " 3 

9m — 



«n(Ql - Qn+l)" 1 - 2 
n-1 



(Ql - Qn+l) k 

, m — 3, 4, . . . , n 



for u n = Qi - 3Q n + ^2 (Qk+i - 3Qfc) 
Proof. Let 



Q„-2 



fc=2 



i — k 



1 

1 











1 , <(Q„-3Q„_i) _ Q„_i <(Q„-i-3Q„- 2 ) _ Qn-2 

2 "/• 



2u„ 2 
Q„-3Q„-i 

"re 
1 



2Mj! 

Qn_-3Q„ 

"re 








<(Q 3 -3Q 2 ) _ % 

2u n 2 
Q3-3Q2 
u n 






1 



n-1 



and G = diag(2, u n ) where u n = Qi - 3Q„ + (<9fe+i - 3Q fe ) ( q~^~ ) 

fc=2 

n n — /c 

and = ^ Q fc ( q~^~7 ^ • Then we obtain 



n — k 



fc=2 



isTQLS = G © A 



where G © A is the direct sum of the matrices G and A. If W = LS, then we 
have 

Q- 1 = WiG' 1 ®A- X )K. 
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Since the matrix <Q> is circulant, the inverse matrix Q 1 is circulant from 
Lemma 1.1 [1, p. 9791]. Let 

QT 1 = circ(q 1 ,q 2 ,...,q n )- 
Since the last row of the matrix W is 

„ , On — 3Q n -i Qn-1 — 3Q n ^2 Qn-2 — 3Q n -3 Q4 — 3Qz Q3 — 3(^2 

) 5 ; ■ * ■ ; 1 

li n ii n ii n li n W n 

the last row elements of the matrix Q _1 are 



<?2 = ( -3 + ^(Q„-fe+l - 3Q„_fc) ( 



fc-i 



93 = - 



fe=i 

q 3 - 3Q 2 



(Qi-Q„+i) fe 



u„(Qi - Q 



n+l, 



2(Q 3 - 3Q 2 ) J_ A (Q 5 -fc - 3Q 4 -fe)(Q n - 2)* 



«n(Ql - Qn+l) Un (Ql - Qn+l) k 



1 ( — 3Q 2 (Qe-k — 3Q 5 _fc)(Q„ — 2)* 

95 ~ ro,_o_.^"2^- 



«n \ (Ql - Qn+l) ^ (Ql - Qn+l) fc 



2^ {Qb-k — 3Q4-fe)(Qn — 2) A 

" ^ (Qi-Q„+i) fe 



1 [ (Q_" -fc+l — 3Q„_fc)(Q„ — 2) A 



(Ql-Qn+l) fc 



2 {Qn-k 3Q„_fe_i)(Q„ 2) fe 1 



fe 



, , (Ql - Qn+l) 

(Qn- fc-l — 3Q„_fe_2)(Qn — 2) fe ~ 1 



fc=l 



(Qi-Q n+ i) fc 



1 / 1 8(Q„ - 2)»" d , (Q 

n — k-\-2 

— 3Q„_fc+l)(Qn — 2) 



u„ ^ (Qi - Q„+i)"- 2 ^ (Qi - Q„+i) fc 
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where u n = Qi-3Q„+^ (Qk+i - 3Q k ) ( q 1 < ^"q" ii+1 ) • Now, let us rearrange 

k=2 

the qiS (i > 4). Hence we obtain 

1 /2Q3-6Q2-Q4 + 3Q3 (Q 3 -3Q 2 )(Q„-2y 
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\ Ql ~ Qn+1 

( 



(Ql-Qn+l) 2 

\ 

-Q 4 + 2Q 3 + Q 2 + 3Qi - Q 2 -4(Q n -2) 



Ql — Qn+1 



(Qi - Q„+i) 2 



V 

4(Q„ - 2) 
U n (Ql - Qn+lY 



J 



Q5 = — 



1 



Ql — Qn+1 
Qn 2 



((-Q 5 + 2Q 4 + Qa) + 3(Q 4 - 2Q 3 - Qa)) 



(Ql — Qn+1. 

4(Q„ ~ 2) 2 
(Qi - Qn+1) 3 



5((-Q 4 + 2Q 3 + Qa) + 3Qi - Q 2 ) + ^3 











If we formulate general case i > 3, then we have 

4(Q„-2)"- 1 



(Ql-Qn+l)' 



Thus 



j_ . / 8(g w -2)"~ 3 

Un C * rC l (Qi - Qn+i)"" 2 



n— 3 

E 

fc=i 



(Qn-fc+2 — 3Qn-k+l)(Qn — 2) 



fc-1 



(Ql-Qn+l)* 



n-2 



-3 + ^(Qn-fc+i — 3Q„_fe) 



/ (Qn - 2) 



fe-1 



fc=l 



\(Ql - Qn+l) k J ' Ql - Qn 



i+l 



4(Q„ - 2) 4(Q n - 2) 2 4(Q„ - 2) 



n-3 



(Qi - Q„+i) 2 ' (Qi - Q„+i) 3 ' ■ ' ■ ' (Qi - Q„ + i)"- 2 
by the nice properties fo Pell-Lucas numbers. ■ 



Corollary 10 Let the matrix K be as in ^Tty. Then its inverse is block matrix 

l= 1 
D E 
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where the matrices 



n _ I Qn Qn-l Qn-2 Q2 \ 

^ \ 2 2 2 2 / 



T 



2 2 2 •••2 /(n-l)xl 

cmd S is an (n— 1) x (n— 1) Hankel matrix that it has the first row [P n -i, P n -2, ■ ■ ■ , Pi] 
and the last column [Pi , 0, . . . , 0] . 



Proof. The matrix K 1 is obtained easily by applying elementary row 
operations to the augmented matrix [K : I n }. m 
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